Non-Fermi liquid and pairing in electron-doped cuprates 
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We study the normal state and pairing instability in electron-doped cuprates near optimal doping. 
We show that the fermionic self-energy has a non-Fermi liquid form leading to peculiar frequency 
dependencies of the conductivity and the Raman response. We solve the pairing problem and 
demonstrate that Tc is determined by the curvature of the Fermi surface, and the pairing gap A(k, co) 
is strongly non-monotonic along the Fermi surface. The normal state frequency dependencies, the 
value of Tc ~ lOK and the k— dependence of the gap agree with the experiment. 
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Physics of high-temperature superconductors contin- 
ues to attract considerable attention. The past research 
primarily focused on hole-doped cuprates, where a com- 

' plex phase diagram was found with antiferromagnetic 
and superconducting phases separated by a pseudogap 

\ region, and where the normal (non-pseudogap) phase dis- 
plays a prominent non-Fermi-liquid behavior. 

] Recently there has been a surge of interest in electron- 
doped cuprates, Nd2-;rCea;Cu04 and Pr2_a;Cea;Cu04. 

[ Their phase diagram |l| has a wider region of anti- 
ferromagnetism (AFM); superconducting Tc{x) forms a 

[ dome above the antiferromagnetic quantum-critical point 
(QCP) [2|, much Hke in heavy fermion materials [3]; and 

] reported pseudogap behavior Q tracks the Neel tempera- 
ture, and so is Hkely just a property of a quasi-2D antifer- 

[ romagnet [5;, 0| . Optical data show [^] the sharp charge- 
transfer gap near 2eV at small x melting away near the 
optimal doping x ~ 0.15, where the superconductivity 
appears. We infer from these data and the experimental 

[ phase diagram that near optimal doping the Mott physics 
is not crucial, and the system behavior can be grasped 
by considering onlylow-energy electrons that interact via 
collective bosons [31 • Electron-boson models have been 

I previously appHed to hole-doped cuprates 0, 0|, where 
Mott physics is more pronounced. They must there- 
fore be tested on electron-doped materials. In particular, 
they should explain a much smaller Tc in electron-doped 
cuprates despite almost the same strength of the Hub- 
bard interaction as in hole-dopes materials 

I Confinement of the superconductivity to the vicinity of 
an antiferromagnetic QCP implies that the most Hkely 
candidates for the pairing bosons are the spin fluctua- 
tions with momenta near the antiferromagnetic vector 

' Q = (tt, tt). RPA studies suggest ^ that AFM appears 
close to the doping when 2kF coincides with Q, i.e. when 
the Fermi surface (FS) touches the boundary of the AFM 
Brillouin zone at ki? = (7r/2,7r/2) (see Fig. Near this 
point the free-fermion static susceptibility x(Q) rapidly 
increases, leading to the Stoner instability. 

Below we set Q — 2kF and consider a spin-fluctuation 
model for the 2kF instability that involves fermions near 
ki? and their collective spin fluctuations near Q. Such a 



model qualitatively differs from the spin-fermion model 
for hole-doped cuprates 0|, where relevant fermionic 
states (hot spots) are located far from the zone diago- 
nals. In that case, the spin-fluctuation exchange led to 
the d2.2_j^2-wave pairing [21 with the gap maxima at the 
hot spots When the dominant interaction is between 
the nodal fermions, it is a priori unclear to what degree 
the c?j,2_j,2-wave pairing survives. 

In this paper we address this and other issues for 2k p 
instability. We conflrm earher result J^] that the inter- 
action with low-energy spin fluctuations destroys Fermi 
liquid behavior for diagonal fermions. We show that 
this causes a peculiar behavior of the optical conduc- 
tivity (j{ll>) over a wide frequency range, and 
a nearly frequency-independent Raman response in the 
B2g channel. Both of the results agree with the experi- 
ment 0, 0| . We then use the normal state results as an 
input and consider for the first time the spin-mediated 
pairing near 2kF instability. We demonstrate that Tc is 
finite at the QCP and the pairing is dominated by non 
Fermi-Hquid fermions. We show that the pairing prob- 
lem is very peculiar as now the curvature of the fermionic 
dispersion plays a major role. We found Tc ~ lOK for 
the same parameters as for hole-doped materials. We 
also found non-monotonic variation of the d— wave gap 
A{k,uj) along the FS (see FigEJ. Non-monotonic vari- 
ation of the gap agrees with the numerical studies 
and with the earher BCS-type study but our so- 
lution is very different from the BCS one. The non- 
monotonic d— wave gap has been detected experimentally 
at smaller x and has been traced to the location of the 
hot spots 0, 0|. Our results show that near a mag- 
netic QCP, the positions of the gap maxima deviate from 
the hot spots, and remain fixed even when the hot spots 
merge and then disappear. 

We now provide the details of the calculations. A spin 
fermion model describes low-energy fermions Ck interact- 
ing with their collective bosonic fiuctuations in the spin 
channel Sq: Hint = gJ^w.q'^^+'i.a'^apci^pSq (for internal 
consistency of the model it is required that g < Ef JJJ)- 
The spin subsystem is described by the static spin suscep- 
tibility which, as before, we assume to have the Ornstein- 
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Figure 1: Left: Fermi surface at the antiferromagnetic QCP 
with 2kF — (7r,7r). Diamond-shaped dashed lines bound the 
magnetic Brillouin zone. The diagonal points of the FS (nodal 
points of the d^2_y2-wave gap) now become "hot". Right: 
Graphic explanation of the attraction in the d^2_y2 chan- 
nel: parts of the Fermi surface on the same side of the zone 
diagonals are on average closer to the Q separation than the 
parts on the opposite sides, leading to attraction in the d^2_y2 
channel (plus and minus are the signs of the d^2_y2 gap). 

Zernike form Xo(q) = Xo/['?^^ + (q ^ Q)^]- Interacting 
fermions are located near kp and kj? -I- Q, and have al- 
most antiparallel velocities: 

e-i.KVF^k^+P'^^kl, eu+Q^ -VFAk^+f3'^Akl, (1) 

where Ak = k — k^?. Factor (3 accounts for the cur- 
vature of the Fermi line. The strategy is to solve self- 
consistently for the normal state bosonic n(q, Q) and 
fermionic S(k, w) self-energies, and then use the full nor- 
mal state propagators as input for the pairing prob- 
lem. We define the fermionic propagator as G{k,uj)~^ = 
i{uj + 'S{'k, uj)) — Ck and bosonic as x(q, ^)^^ — Xo^i'i) + 
2[n(q,c^)+n(2Q-q,c^)]. 

Normal state Fermionic and bosonic self-energies in 
the normal state have been obtained by Altshuler et 
al. fl^ within Eliashberg theory, and we just quote their 
result: at the 2kF antiferromagnetic QCP and k = kp 
along the diagonal, 

n{Q,n) = ^^^\, E(kj.,c.) = c.y^|c.|3/4, (2) 

where g — g^Xo, and luq — [gP/TrNvp)'^- Both 
self-energies differ from the case when the hot spots 
are far from zone diagonals. E.g., n(q, fi) scales as 
instead of the conventional Landau damping form 
n(q, r2) cx This is a consequence of antiparallel ve- 
locities at ki? and ki? -I- Q (observe that the damping 
rate Il{q,fl) cx does not diverge only because the 
curvature /? ^ 0) . Eq. [2l is strictly speaking valid for 
LU > ujQ. For LU < ujQ expressions for the self-energies 
are more involved |13j, and one has to sum up series of 
logarithms to obtain the proper exponents. Altshuler et 
al. found T,{lu) luqlu^'"', where a ~ 0.14. We con- 
firmed their result. This low-frequency behavior, how- 
ever, is less important for practical purposes, since we 
found that luq ~ lOmeV is quite small 20] . Relevant 
physics, including the pairing, then mostly comes from 
frequencies u > loq, where Eqs. ^ are valid. 




Figure 2: Normal state conductivity as a function of frequency 
shows a scaling behavior for luq < < 40(j;o. Left: |(j(ci;)|, the 
behavior is indistinguishable from the o-{uj) ~ oj""'^'* depen- 
dence. Right: arctan 9(j(aj)/5R(T(tj), an almost constant value 
means that both 9(T(a;) and 5R(T(a;) scale as lo~'^'^'^. 



For the FS momenta away from the diagonal {ky ^ 0, 
see Fig. QJ, the non- Fermi liquid form of the self-energy 
(EJ survives up to fcy^max ~ w^/^, and crosses over to the 
Fermi liquid form at larger ky. At the same time, for k 
transverse to the FS the self-energy is regular already at 
ky = 0: E(fca,, 0) oc kx. 

Optical conductivity We used the result for the self- 
energy to compute complex conductivity <j{iu) from the 
Kubo formula appropriate when the self-energy mainly 
depends on frequency. We verified that the dominant 
contribution to <jiuj) comes from nodal regions where 
S(ti;) has a non- Fermi-liquid form At very small fre- 
quencies, LU < ujQ, S(a;) > cu, and both Scr and Sftcr have 
the same power-law behavior: (y{u!) oc ^ w^*^'^. 

At larger frequencies, S(w) < lu, and one should not 
generally expect 3?cr and 3cr to scale with each other. 
Surprisingly, the scaling behavior extends, with almost 
the same power, a{uj) oc a;~°-^^ up to ~ AOluq (see 
Fig. EJ. Such power-law behavior is not a sign of a 
true scaling, but rather a consequence of the fact that 
'S{uj)/lu oc (luq/luY^^ is a slow decaying function. The 
Lu'f behavior of conductivity with 7 w 0.68 has been ob- 
served in Pri.85Ceo.i5Cu04 below 400meV Both 
the exponent and the experimental frequency range are 
quite consistent with our results. Note that a very simi- 
lar behavior of the conductivity at intermediate energies 
has been earlier observed in hole-doped cuprates |23 |. 

Raman response We also considered the Raman scat- 
tering in the B2g channel. Raman intensity is propor- 
tional to the imaginary part of the polarization bubble 
weighted with the geometry dependent Raman form fac- 
tors 7fe. In the B2g geometry, 7(k) oc sin fca; sin fcy , such 
that the signal comes primarily from around the diago- 
nal direction. Unlike the current vertex for conductivity, 
the Raman vertex is renormalized by the interaction even 
when S(a;) depends only on frequency. Because 7(k) is a 
scalar and doesn't change sign under k ^ k-t- Q, we can 



Figure 3: Pairing instability temperature in units of 
the coupling constant g vs. the curvature parameter r — 
g[3^/-KVp. Solid line gives the result of the numeric solu- 
tion, which requires increasingly long computational time 
as r — > 0. Dashed line is the analytic small-r form Tc = 
0.955r/(l + 5r2/5)^ 

approximate it by a constant jCkp). The renormaliza- 
tion of the Raman vertex then coincides with that of the 
density vertex, and is related to the self-energy by the 
Ward identity: the full 7fuii(k,cj) = 7(k)(l + d^Y.{u;)). 
Evaluating the Raman bubble and substituting 7fuii into 
it, we then obtain 

Rb,, (uj) - k„{^) (1 + d^n^)f /I (3) 

-j- Zj yui j 

Substituting the forms of the self-energy and A:j^,max('^) ~ 
u)" for Lu < loq, fcj,.inax(^^) ~ w^^^ for LO > loq, we find that 
at r = 0, RB2g{^) is flat: it is a constant at w < loq, 
and very slowly crosses over to the (cj/wq)"^^'' behav- 
ior at a; > Wo- The flat form of the Raman intensity at 
frequencies uj > ujq is consistent with the experimental 
data The experimental behavior at w < is dom- 
inated by temperature effects which we do not consider. 

Pairing instability We now apply normal state results 
to study the pairing instability near a magnetic QCP. As 
we said in the introduction, it is a priori unclear whether 
d— wave pairing survives when the hot spots merge, as the 
strongest pairing interaction involves quasiparticles for 
which the dx2-y2 superconducting gap vanishes. How- 
ever, separation between near-nodal fermions with oppo- 
site signs of the d^2_y2 gap is on average closer to Q than 
that between fermions with the same sign of the gap (see 
Fig. n]b). So the spin-fluctuation exchange still leads to 
attraction in the d— wave channel. Furthermore, as the 
gap equation relates A(k) near and k^ + Q, where 
both gaps are linear in the deviation from the nodal di- 
rection, the gap equation becomes an equation on the 
slope of A(k) and does not contain any extra smallness. 

The linearized gap equation at QCP is obtained using 
Eliashberg technique for collective-mode mediated pair- 
ing We verified that the dominant contribution to 
the pairing comes from uj > ujq, where T,{uj) can be ap- 
proximated by the uj^^^ form. To single out the role of the 



Figure 4: The pairing vertex ^{ky,Lij„) for r = 0.05 vs 
cj„/(27rTc) and ky/qo {qo is defined in the text). Observe 
that ^{ky,uj„) is non-monotonic in ky. 

curvature, we rewrite loq as loq — ujir, where uji — g/(47r) 
accounts solely for the interaction, and the dimension- 
less r = g(3'^/TTVp ~ g/Ep parameterizes the effect 
of the curvature. Introducing dimensionless variables 
k/qo, uj/uji, where go = g/ {■nvpr'^/^) ^ kFig/Ep)^^*, 
and assuming that the gap has d^2_y2 symmetry, i.e., 
A(kp) — — A(ki? -I- Q), we obtain an integral equation 
for the pairing vertex $(fcj^,w„) near a single hot spot: 

r ^{k'y,U^)dk'y 

J iky-k'y)^+ry^ki + k'y^)^ + ./\^j;:^\' 

The momentum dependence in the kernel in Q comes 
from (k — k')^ after the substitution kx = fi'^ky/vp and 
k'^ = —(i'^k'^jvp. The pairing vertex is related to A as 
<i> — Aiojiuj -I- S(ti^)). A (i2,2_j^2-wave solution of I^J is 
antisymmetric in ky and symmetric in w |2^ . 

At vanishing r, one can neglect the curvature of the 
fermionic dispersion in the gap equation (the r^^"^ term 
in the denominator in Q). Without this term, the kernel 
in Q is non-Fredholmian, and its antisymmetric-in-fc so- 
lution is a monotonic $(fcj,,Ci;„) = fcj,$((jj„). Substituting 
this form into Q, we obtain for <i>(a;„): 

$(w„) = \.T y: ^i^'^M — ^ 

m^n |t^„-m|^/'* f|Wm| -|- Wq' | 1 3/4 J 

(5) 

Expectedly (since we neglected the curvature) , this equa- 
tion is the same as in the Eliashberg theory for the pair- 
ing, mediated by a bosonic susceptibility x(a;) oc cj^^/** 
[2J|. Numerical solution of (UJ gives Td ~ 6ujo = 6wir, 
i.e. Tc ~ r. We emphasize that the r ^ limit is not 
weak coupling (self-energy cannot be neglected) despite 
the factor r^/^ in the numerator of 

When the curvature of the fermionic dispersion is in- 
cluded, Eq. (21 becomes Fredholmian, and a numer- 
ical solution can be obtained using the Nystrom dis- 
cretization method [23. We found that the effect of 
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the curvature on Tc is very strong: above r > 0.001, 
the actual Tc rapidly becomes much smaller than Td. 
For r ~ 0.1, - O.OlT^i. We plotted Tc(r) in Fig. 
El We see that over a wide range 0.01 < r < 0.1, 
Tc ~ O.OOGwi K, 0.0005^ forms a plateau and weakly de- 
pends on r. Using the same g ~ 1.6eV as for hole-doped 
cuprates |2^, we obtained in this range Tc ^ lOK. Using 
further the previous estimate of wq = wir ~ lOmeV we 
obtain r ~ 0.08, which is within the region where Tc is al- 
most a constant. Out Tc ~ lOif at a mag netic QCP is in 
agreement with the experiment [l|, [a, • Note that the 
experimental Tc increases in the antiferromagnetic phase 
up to 20K, and only then drops at smaller dopings. 

The huge discrepancy between Tc and Td can also be 
understood analytically, by expanding Tc in r beyond 
the 0{r) term. We found that the expansion is rather 
non-trivial and holds in fractional powers of r2/5: Tc « 
Tci/ (1 + ar^/^)^. The prefactor a > cannot be obtained 
analytically, but the fit to the numerical data gives a ~ 5. 
As (1 -I- 5r^/^)* rapidly increases with r (it is 80 for 
r = 0.1), this explains why Tc is much smaller than Td. 

We also found that at finite r the gap is a non- 
monotonic function of ky-. it is linear at small fcy, passes 
through a maximum at a fco and then falls off (see Fig. 
^ . Position of the maximum depends on r and appears 
uncorrelated with the positions of hot spots, which at 
the QCP are on the zone diagonals. At small r we found 
analytically that ko{r) ~ kpr^^^ ~ kp{g/ Ep)^/^ . As r is 
small, fco < kp and the gap is confined to the vicinity of 
the zone diagonals. 

At deviations from the QCP towards larger dopings, 
into a paramagnetic phase, Tc decreases and eventually 
disappears. The value of fco, however, does not track the 
decrease of Tc, i.e., the d— wave gap extends over a finite 
momentum range along the FS even in the overdoped ma- 
terials. At deviations into the antiferromagnetic phase, 
the FS evolves into hole and electron pockets, and the 
locations of fco gradually approach the hot spots. 

Our results for Tc and the gap survive even when the 
magnetic correlation length ^ remains finite at the doping 
where 2kp — Q, i.e., when antiferromagnetic QCP is 
shifted to lower dopings, when the hot spots are away 
from zone diagonals. We found that the modifications of 
our results are small as long as Tc > J^~^ where J is the 
exchange interaction. 

Summary We considered the normal state proper- 
ties and pairing near a 2kp antiferromagnetic QCP and 
applied the results to electron-doped cuprates at optimal 
doping. We found that the breakdown of the Fermi-liquid 
description at QCP leads to peculiar frequency depen- 
dences of the conductivity and the i?2g Raman response. 
We found that the superconducting Tc remains finite at 
the QCP and is of order lOK for the same spin-fermion 
coupling as for hole-doped cuprates. The pairing gap at 
QCP has da;2_j,2 symmetry, but is highly anisotropic and 
confined to momenta near the zone diagonals. 
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